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Abstract 

The nature and symmetry of the superconducting gap function in the noncentrosymmetric super- 
conductor (NCS) Li2Pt^B, even many years after its discovery, appears to be full of contradictions. 
In this letter based on the existing band structure calculations we find that owing to the consider- 
able nesting near the Fermi surface and the enhanced d-character of the relevant bands that cross 
the Fermi level, the system gets somewhat strongly correlated. Considering the effect of the onsite 
Coulomb repulsion on the pairing potential perturbatively, we extract possible superconducting 
transition. The strong normal spin fluctuation gives rise to a singlet dominant gap function with 
accompanying sign change. Thus our theory predicts a s± wave gap function with line nodes as 
the most promising candidate in the superconducting state. 

PACS numbers: 74.20.Mn,74.20.Rp,74.70.-b, 74.90.+n 
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The occurrence of superconductivity in compounds without spatial inversion symmetry 
is one of the most active fields of research now a days. Inversion symmetry breaking leads to 
many new interesting effects in the superconducting state. The discovery of superconduc- 
tivity in Li2Pd^B [1] and subsequently the experiments in the pseudo-binary complete solid 
solution Li2{Pdi_xPix)zB , a; = ~ 1 [2] attracted much attention. Many experimental as 
well as theoretical works reported since then. The end compounds Li2Pd^B [x — 0) and 
Li2PtzB [x — 1) [3-8] were also studied intensively and compared. It is now established that 
the superconductivity in Li2Pd^B is phonon mediated s-wave type. The presence of Hebel- 
Slichter peak in the Nuclear spin-lattice relaxation rate measurement [3, 4], low temperature 
behavior of the specific heat [5, 6], penetration depth [8] etc. as well as NMR Knight shift 
data strongly support this conclusion. On the other hand, the nature and symmetry of the 
gap function of the compound Li2PHB is still debatable. Similar experiments performed on 
this compound [4, 6, 8] suggest the presence of line nodes in the superconducting state. The 
NMR Knight shift, often used to distinguish the spin state of superconductivity between 
singlet and triplet, is almost temperature independent even below Tc. This behavior is also 
interesting and deserve special attention. 

In this letter, based on the existing band structure calculations [9, 10] we find that there 
exists considerable nesting between the Fermi surfaces and enhanced d-character of the 
relevant bands that cross the Fermi level. These two effects lead to a stronger electron 
correlation in Li2Pt^B than in Li2PdzB [11]. By treating this correlation perturbatively we 
estimate the effect of spin fluctuations on the stability of superconducting state [12]. Owing 
to the large usual spin fluctuation the singlet gap function becomes much stronger than the 
triplet gap function. The singlet gap function belongs to Ai representation with sign change 
between two branches of the Fermi surface. Thus our theory suggests a singlet dominant s± 
wave gap function with line nodes as the most promising candidate for the superconducting 
state of Li2Pt'iB. This prediction also explains most of the experiments. We also calculate 
the behavior of uniform spin susceptibility below Tc and comment on the apparent mismatch 
with the experiment. 

The crystal structure of the compound Li2Pt^B is simple cubic (with point group O) 
and isostructural with the compound Li2Pd^B. The only difference between them is in the 
mass of the central Pt and Pd atoms. However this gives rise to some significant observable 
effects [9, 10]. For Pt compound there is an enhancement of d-character of the bands that 
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FIG. 1: (Color online) Signature of nesting. (a) Momentum dependence of interband and intraband 
susceptibilities showing nesting at Q = (vr, vr, 7r).(b)The comparison of the relative magnitudes of 
normal and anomalous spin fluctuations. The charge fluctuation Xcc and normal spin fluctuation 
\{Xxx + Xyy + Xzz) contribute strongly compared to the anomalous spin fluctuations. The spin 
susceptibility also get enhanced at the R point (details in text). 

cross the Fermi level. This enhancement of the d-character is reflected in the increased 
DOS at the Fermi level. Considering these, we can construct a minimal model Hamiltonian 
{H = Hq + Hi) of Li2Pt'iB which is given by the Hubbard model with an antisymmetric 
spin-orbit (SO) coupling term, where 

^0 = XI ([^k - lA + gk-0-)^^, 4aCka', (1) 
\icrcT' 

and Hi = U J2i ''^it^ii- Here Cko- and cj^^ denotes the annihilation and creation operators 
of an electron with momentum k and spin a. is the dispersion of electrons and ^ the 
chemical potential, gk = — g-k denotes the effective anti-symmetric SO coupling which 
breaks the inversion symmetry. In i^i, U is the screened on-site interaction. The dispersion 
of electrons is constructed by the tight-binding method including upto fourth-neighbor 
hopping in the three-dimensional simple cubic lattice. 

= 2ti{cos{kx) + cos{ky) + cos{kz)) + 4:t2{cos{kx) cos{ky) + cos{ky) cos^k^) 
+ cos{kz) cos{kx)) + 8^3 cos{krc) cos{ky) cos{kz) + 2t4{cos{2kx) + cos{2ky) + cos(2A;2))(2) 

The SO coupling term appropriate for the point group is given as gk = 
g{sm{kx),sm{ky),sm{kz)). The values of the parameters (ti,t2,t3,t4, g, ^) are chosen to be 
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FIG. 2: (Color online) The variation of maximum eigenvalue with U. 

(1.0, —0.03, —0.88, —0.03, 0.5, 0.02) as Fermi surface obtained by the band structure calcula- 
tion are reproduced. One can diagonalize Hq to get the eigen-energies ek± = ± |gk| — A*- 
The spin degeneracy is now removed and ek± gives us the energy of the two helicity bands. 
The Fermi surfaces corresponding to the helically splitted bands consist of three major 
branches. One electron pocket around the F point, one hole pocket around the R point 
and other remaining parts forming a cage like structure with neck and mouth along T — X 
direction. The corresponding parts of the Fermi surfaces of the different helicity bands are 
shifted from each other depending on the magnitude of g. Owing to the smallness of the 
parameters ^2,^4 and fi there appears a large nesting with the nesting vector Q = (7r,7r,7r) 
connecting between the cage like larger portion of the Fermi surface of the negative helicity 
band ^k- and the similar cage like Fermi surface of the positive helicity band £:k+- There- 
fore this nesting gives rise to a sharp pick at the R point in the momentum dependence 
of x+~ (q) = 8^ '^'''^£1^+ /-e'l^''^^ where nesting condition is satisfied. There also ex- 
ists partial nesting between Fermi surfaces around F— and R-points of both helicity bands 
but for positive helicity band it's larger as shown in Fig. 1(a). The most general form 
of the matrix gap function is Ak = [\l'(k)(To + d(k).o"] i(3"y. Here ^'(k) is the singlet gap 
function and d(k) is the triplet d-vector. In NCS superconductor triplet component with 
|d(k).gk| = |d(k)||gk| survive the pinning from the antisymmetric SO coupling [13]. So one 
can write d(k) = 0(k)gk where 0(k) having the same symmetry of momentum dependence 
as ^^(k). With all these we can define the normal (j(k, iun) and anomalous F{k, iujn) matrix 
Green's functions as below, 

^(k,^^^) = G+(k,iu;„)(3-o+G'_(k,iw„)gk.a, F(k,iu;„) = [F+(k, iw„)6-o + F_(k, iw„)gk.(3-] io-^^ 

(3) 
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here gk = gk/|gk|- G± and F± are given as, 
G±(k,za;„) = - — — — 3— ± — — — 3— , F±(k, ' - 



here Ak± = ^(k) ± 0(k)|gk| and £'k± = y^£^_|_ + A^_,_. Within the weak couphng theory of 
superconductivity only static susceptibihty is required. We start by defining the dynamical 
susceptibility as, 

l-l/T 

^0 

here (...) denotes thermal average, 7^ imaginary time ordering and Q„ are the Bosonic 
Matsubara frequencies. The charge (spin) operators with wave vector q is defined as, 

k(T k(T<j' 

With all these, the matrix elements of the static spin susceptibilities Xa/9(q) for a,/? — 
c, X, y, z is found to be, 

where f{e) is the Fermi distribution function and the function F^^ is obtained as, 

(k; q) = 6a,l3il - CCgk-gk+q) + CC(^ka^k+q/3 + ^k/3^k+qa) " e„;3^i({^k+q7 - C^k7)- (7) 

Similarly the charge fluctuation in the normal state i.e. Xcc(q) is obtained with the re- 
placement of F^||'(k;q) by r^^(k;q) where r^^(k;q) = 1 + '^Cgk-gk+q- The susceptibilities 
between spin and charge operators Xca(q) and Xac(q) all vanishes for the static case. We 
calculate all the susceptibility components and examine the property of spin fluctuations. 
The usual spin fluctuation ^{xxx + Xyy + Xzz) with the momentum dependence of q^— type 
is also present in the centrosymmetric cubic system. Other symmetric spin fluctuations 

zz ~X,xx X.yyjj \X,xx 

— Xyy) and {Xai3 + Xi3a) with a ^ P having momentum dependence 
~ Qx ~ (lyj ^x ~ (ly and qaqiii^a 7^ /?)— types respectively arc special to the cubic noncen- 
trosymmetric case [12]. Along with these the anti-symmetric spin fluctuations i{Xai3 ~ Xi3a) 
with a ^ (3 with momentum dependence q^{j ^ a ^ (3)— type are also present. In Fig. 1(b) 
we compare relative strengths of the charge and usual spin fluctuation together with the 
anomalous spin fluctuations along symmetrical lines. Later we will see that the largeness 
of the usual spin fluctuation is responsible for the largeness of the singlet gap function as 



FIG. 3: (Color online) The singlet gap function at three values; (a) = 0.0, (b)^^ = O.Tvr, and 
{c)kz = vr. Gap function varies from positive- maximum (Red) and to the negative- maximum (blue) 
following the scale attached with each figure. The blue line with + sign shows the Fermi surface of 
the positive helicity band while the dotted blue line denotes the nodes of the gap function (details 
in text). 

triplet gap function whose magnitude is much smaller than the singlet one are induced by 
the antisymmetric spin fluctuations. 



Treating the interaction term perturbatively and following the standard procedure [14-16] 
we arrive at the following superconducting gap equation, 



-E 



(8) 



f VUq) Kx(q) K,(q) K.(q) \ 
V^xs(q) Vxx{q) Vxy{q) V^ziq) 
Vysiq) Vy^iq) Vyyiq) Vy,{q) 

V vUq) Kx(q) v.yiq) v^q) J 

here V^^ with {(, t]) = [s, x, y, z) denotes the pairing potential arising from the corresponding 
fluctuation exchange and they are expressed as below. 



( ^(k) \ 
4(k) 

dyiX) 



( J-.(k-q)\ 
^x(k-q) 

^y{^ - q) 

V J-.(k-q)/ 



Ks(q) = -u -u^ [xxx(q) + Xyy{.q) + Xz-M) - Xcc(q)] 
%(q) = u"^ [xcc(q) + Xvn^q) + xsi^q) - xcc(q)] 
^.(q) = K,c(q) = -u^ [xM + Xr,c(q)] 

Kc(q) = -'^C«(q) = [Xrisi.'^) - X5r,(q)] , 



(9) 



where C 7^ ^ 7^ ^- Here J^g and J-'q are the contributions from the anomalous Green's 
functions after frequency summation [14]. As we mentioned above the singlet component of 
the gap function arising from usual spin fluctuation dominates over triplet gap function which 
is induced by the small antisymmetric spin fluctuations. Eq. (8) reduces to the eigenvalue 



problem if we work at the transition temperature Tc. We fix % arbitrarily at 0.02ii and solve 
Eq. (8) for maximum eigenvalue. Fig. 2 gives us the critical value of the onsite interaction 
U for superconductivity i.e. Uc — 1.675ii , when maximum eigenvalue becomes unity. 
We thus get the momentum dependence of both singlet and triplet gap functions as the 
eigenfunctions of the maximum eigenvalue. From the momentum dependence of the gap 
function we conclude that the superconductivity belongs to the Ai representation of the 
point group O. 

In Fig. 3 we present the contour plot of the singlet gap function at three kz values in the 
1st Brillouin zone. The singlet gap function changes sign from positive (red) to negative 
(blue) gradually moving from F— to R-point and vanishes completely somewhere in between 
forming the nodal surface. Here we would like to mention that the gap function is strongest 
at either F— or R-points although the nesting is rather weak here. On the other hand the 
cage like portion of the Fermi surface where we have most strongest nesting gives rise to weak 
gap function. It can be understood from a careful observation of Eq. (9). The summation of 
prefactor F"" (k, Q) of susceptibility Xaa{Q) vanishes for the spin-singlet pairing potential 
Vss(Q) at the nesting vector Q = {n,n,n). Because of this, even strong interband nesting 
does not play any role in opening up the gap function. However, the pairing potential forms 
the gap on the Fermi surfaces around F— and R-points, connected by sub-dominant nesting 
of x++(Q). Thus the singlet gap function with opposite signs between these points opens 
up and line nodes can exist in between. In this figure we show the corresponding zeros of 
the positive hehcity Fermi surface by the -|- sign fine and the dotted line denotes the exact 
location where the gap function vanish. In Fig. 3(a) around F— point the gap is positive 
maximum. In Fig. 3(b) the strength of the negative gap function increases and finally, in 
Fig. 3(c) we encounter the maximum negative value of the gap function at the corner R- 
point. Thus the gap function appears to be singlet s± type with accidental line nodes which 
is not allowed by symmetry, rather depend on the three-dimensional geometry of the Fermi 
surface. 

We also calculate the temperature dependence of the susceptibility in the superconducting 
state. Within the weak coupling approximation neglecting the feedback effect we assume 
that the order parameter below follows the BCS temperature dependence A(k, T) = 
A(k, 0) tanh(1.74-i — 1) at every momentum point. Using this gap function we calculate 
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FIG. 4: (Color online) The variation of normalized spin susceptibility (solid red line) with tem- 
perature in the superconducting state. Also shown in the same figure the contributions from 
the Van-Vleck term (dashed magenta line) and the Pauli term of positive helicity band (green 
square)and negative helicity band (blue square), all scaled by Xs(^c)- 

the uniform susceptibility below Tc as follows [17]. 

Xs{T) = [Xo(k, T) + x+(k, T) + X-(k, T)] , (10) 

k 

where xo(k, T) , x_|_(k, T) and x_(k, T) are given as, 

Xo(k.T) ^ 15:1(1 -f^^|-±it±£t) X (!?2!^(M±||^(i)j]. (U) 

X-±(k,r) = l/(12Tcosh'(Bk±/2T)), (12) 

We plot in Fig. 4 the contribution from the temperature independent Val-Vleck term xo{T), 
temperature dependent Pauli terms x+iT),X-(T) and the susceptibihty Xs (T) all normal- 
ized by Xs{Tc) (red line). Although Xs{T)/Xs{Tc) shows excellent agreement with earlier 
works [18] but apparently contradicts the NMR Knight shift data [4]. To explain the con- 
tradiction with the experiment, one can formulate a multi-orbital theory which captures the 
complicated band structure in more detail. Then the large contributions of the Van-Vleck 
term between t2g— and Cg— orbitals for cubic system is expected and this will further re- 
duce the deviation of the normalized susceptibility from normal state below T^. This involves 
somewhat elaborate calculations and we leave this as a future problem. 

In conclusion, we suggest that, in the noncentrosymmetric superconductor Lz2-Pt3-B con- 
siderable d-character of the bands near the Fermi energy and nesting of the Fermi surfaces 
give rise to weak correlation effect which can be treated perturbatively and this give rise to a 
singlet dominated (with negligible triplet component) s± kind of gap function with acciden- 
tal line nodes arising from the Fermi surface geometry. The three-dimensional geometry of 
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the Fermi surface and the nesting of the Fermi surface also play a crucial role in determining 
the nature of the gap function. We propose that angle-resolved photo emission spectroscopy 
and de Hass-van Alphen effect experiments may shed light on this nesting property of the 
Fermi surface and can be useful to study the properties of the superconducting state as 
well. We also calculate the susceptibility below and emphasize the importance of orbital 
degeneracy of d-electron to explain the experimental data. 

The authors are grateful to Y. Yanase for fruitful stimulating discussions during his visit 
in APCTP. 
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